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Abstract We calculate atom-photon resonances in the Wigner-Weisskopf model, ad- 
mitting two photons and choosing a particular coupling function. We also present a rough 
description of the set of resonances in a model for a three-level atom coupled to a scalar- 
photon field. We give a general picture of matter-field resonances these results fit into. 
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I. Introduction 

Last years there has been a renewed interest for atom-radiation states in non-relativis- 
tic quantum electrodynamics. After the rigorous proof^ of the existence, for smaU couphng 
constants, of resonances (singularities of an analytic continuation of the Hamiltonian re- 
solvent), resonances coming from the naked- atom-Hamiltonian eigenvalues, many studies 
have been concentrated on the study of the Hamiltonian spectrum and more particularly 
on the existence of the fundamental state in various models^'^'^'^'^'^, for arbitrary coupling 
constants. This latter question is a non trivial problem since it has long been known that 
negative Hamiltonian eigenvalues may appear when the coupling constant increases^'^. In 
preceding works (see Ref. 10 and references therein), for our part, we were interested in 
a method for calculating resonances for arbitrary values of the coupling constant. It is a 
way of catching the above mentioned particular eigenvalues. The result we obtained can 
be generally stated as follows. The coupling of a discrctc-lcvcl system to a zero-mass 
field does not only shift the level energies into the complex plane. Certainly these energies 
become resonances. But the coupling also creates other poles of the resolvent (or of its 
continuations), which have to be placed on the same footing as the preceding ones. For 
certain values of the parameters of S, or certain values of the {S+ field) coupling constant, 
these latter poles may be eigenvalues of the coupled-system Hamiltonian. (In the paper, 
the word resonance will refer to such a pole or, by extension, to an eigenvalue, when there 
is an obvious possible continuous transition from the one to the other.) 

The presence of these poles is well-known for an atom or a molecule in an envi- 
ronment in which the only emitted or absorbed photons are monochromatic (lasers or 
cavities). Indeed, for a two-level atom with energies cq and ei, and photons with en- 
ergy ei — eo, the upper level is split into two levels by the coupling, when at most one 
photon is considered^ (vacuum-field Rabi splitting). The coupling increases the number 
of Hamiltonian eigenvalues because it splits the degeneracy of each of the eigenvalues 
ei, 2ei — eo, • • • , ei + n(ei — eo), • • • of the uncoupled atom-photon system. This phe- 
nomenon occurs in different but analogous situations: for instance the coupling of an 
exciton to the mode of a cavity^^'^^, and it can be detected by spectroscopic means. It 
is also present in electron-phonon interactions^^'^^. Let us note that the final number of 
atom-photon states is simply a consequence of the combination of discrete atom states 
with discrete photon-states. 

Many papers have studied the coupling of a two-level system to another system which 
has either discrete levels or a continuum of levels, often focusing on the continuous transi- 
tion from the vacuum Rabi splitting to the Fermi's golden rule^^'^^. In the present work, 
a splitting is exhibited also in the continuum case. The eigenstates of the naked atom we 
are used to are thus but a small part of all the possible states. 

In fact, the appearance of new resonances when the system S is coupled to a field is 
a general phenomenon whose explanation is given in Section II. It is, of course, not due to 
the smallness of the photon-state width. A numerical method, explained and illustrated in 
previous publications, enabled us to study and calculate these resonances in various simple 
models: S was either a harmonic oscillator^^ or a two-level atom. But in this latter case, 
we almost always limited ourselves to considering physical states with at most one photon. 

A reason to carry on with the study is that the new resonances appearing with the 
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coupling might play a role for systems having a large spatial extension: for instance Ryd- 
berg atoms or, a case maybe more important, molecular orbitals of big molecules. Indeed, 
in a very simplified model for an interaction of S with the field, a model in which the 
spatial extension S of S may be varied {S is a charged harmonic oscillator whose mass 
and spring constant may be varied), we noted that some among the poles we speak about 
come to the negative real axis when S gets large enough^^. These poles then correspond 
to stable states and are therefore important states. Their wave function can be written: 
it is a mixing of electron and photon states. Analogous phenomena may be expected for 
more complex extended systems. 

In this perspective, it is important to be able to start the study of the resonances in 
the two following situations: a system with more than two levels and, in the case of two 
levels, a system with several (non monochromatic) photons. Section III deals with this 
second question. To the best of our knowledge, it is tackled in the literature only when the 
fundamental state question is discussed'*. Here we calculate some mixed (or hybrid) states 
with several photons, or the resonances which correspond to them. Section IV tackles the 
same problem for three-level atoms, but since the situation is more complicated, we limit 
ourselves to a qualitative description of the numerous resonances. Both studies lead to the 
reasonable conjecture that the number of resonances should be the product of the number 
of atom states by the number of independent radiation states actually coupled to the atom. 

II. Notation of the resonances 

We need a precise notation for the different poles. 

It will follow from a very general argument, which also explains why every atomic level 
should give rise to a double infinity of resonances. Let us consider a material quantum 

system 5' whose Hamiltonian, Hs, has eigenstatcs | 0), | 1), | 2),---, with energies 
Co, ei, 62, ■ ■ ■. Let us suppose that this system is coupled to the field of a massless boson, 
here the photon. We denote the state-space of the {S + field)-system by S. 

Let us consider a very general form for the Hamiltonian of the coupled system: 



where i^rad is the energy operator for the photon field and A V represents the coupling of 
S to the field. Let us introduce the auxiliary Hamiltonian. 



where jj is a parameter which may be zero or positive. 

If H — and A = 0, the energy levels cq, ei, ■ ■ ■ are infinitely degenerated (in £), as the 
number of photons accompanying state | m) may be any integer and, moreover, the 
dimension of the space of possible photon states is infinite. Thus the dimension of the 
eigenspace Tii associated with the eigenvalue of H{0, 0) is infinite. Let us underline the 
fact that the degeneracy we speak about here is different from the one we mentioned in 
the second paragraph of the introduction. 

The idea underlying our present study, as well as the preceding ones, is to perturb H{0, 0) 
with respect to A and fi successively. A priori, the perturbation with respect to A removes 



H{X) = i?s 1 + 1 (8) i^rad + A F 



(2.1) 



H{X, IJ,) = Hs^l + fMl^ i?rad + A F 



(2.2) 
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the degeneracy, although it may be only partially, in some particular cases. This leads to 
the following notation. 

Notation 1 (^ii(A), Ci,2(A),--- denote the JN'^ (= IN) eigenvalues of H{X,0) which 
tend to Ci when A tends to 0. 

The choice of the second index, or an equivalent choice, will be described in the 
different contexts. 

The perturbation with respect to /j, then leads us to the following notation for the 
resonances we are interested in: 

Notation 2 When we have selected, in a way which remains to be defined, one pole 
among those of the resolvent of H{X,iJ,) which tend to Ci,jW when /i tends to 0, we denote 
it by Zij{X,iJ,) . 

There are two reasons for having to make a selection. The first one is that if the Ci,j 's 
remain degenerated, and this will be the case in the model of Section III (see Proposition 
3.2), then the degeneracy may be removed for // 7^ 0; in other words, for fi ^ 0, there 
may be several resonances (or eigenvalues) going to the same when /i goes to 0. (See 
a more precise description at the beginning of Section III. 4. 2.) The one we call Zij{X,iJ,) 
will be a particular one, selected in a way which will be made precise in the context. The 
second reason is that the matrix elements of [z — H{X, /u)]~^ are multivalued functions of 
z for fi ^ 0. Therefore, two poles of two different determinations may have the same limit 
when n goes to 0. However, we assume that there is only one pole tending to Qj for a 
given determination, once the first choice has been made. This will enable us to get the 
Zij{X, //)'s by numerical calculus, starting from their germs Ci,i(A)'s. The notations of the 
poles of the different branches of the resolvent matrix-elements will be made precise later 
on, by adding an upper index which refers to the determination (see Section III. 4. 2. 3). 

Ill Hybrid states for a tvi^o-level atom in the Wigner-Weisskopf model 

Two resonances (distinct or not, see below) are well-known in the Wigner-Weisskopf 
model, or in the Friedrichs model, which describe a two-level atom coupled to the field of a 
massless scalar boson^'^'^. These resonances can be seen without considering two- photon 
states. One can be associated, for small coupling constants, to the excited state of the 
atom; the other one is an eigenvalue corresponding to a stable state which differs from 
the unperturbed fundamental state. This eigenvalue appears when the coupling constant 
gets large enough. (Up to recently, it was not clear whether this two resonances were 
two occurrences of the same resonance or not. In fact this depends on the parameters of 
the physical system, but in any case it can be shown that there are actually two different 
resonanccs^*^, for a given value of the coupling constant). In this section we want to study 
other resonances by taking into account several photons. The existence of such states 
is already alluded to in Ref. 4, where Theorem 2.2 shows that the ground state of the 
coupled system must take several-photon states into account. In a particular coupling, we 
will construct two of these states or resonances. It is the subject of Section III. 4. 2. In 
accordance with the general notations introduced in Section II, they will be denoted by 
zo,2{^, 1) and -2i,2(A, 1), tending respectively to the energies and 1 of the naked atom. 
The index 2 will be explained later on. 
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111.1 The model and some notations 

The atom state-space is C^. The fundamental state is | 0) :— (1,0), with energy 

Co = and the excited state is | 1) := (0, 1), with energy ei = 1. In the { | 0), | 1)} basis, 

/O 1 \ °° 
the annihilation operator is a = „ „ . The field state-space is ^ := ® as usual, 

\ U U y n=0 

with J^o ■■= C and J^n ■= L^(IR)^". Let g be in L2(]R) with H^H^ = 1 and c{g), c*{g) 

be the operators annihilating and creating a (scalar) photon in state g. Let i^rad be the 

energy operator in T. 

The Hamiltonian of the model is 

H{\) :=a*a0l + l(g)iyrad + A(a* ® c{g) ^ a ® c* {g)) (3.1) 

In order to have the mixed states appear readily, we introduce 

H{\n) :=o*a(g)l + // I ® H,^^ + \{a* ® c{g) ^ a ® c* {g)) (3.2) 

Note that if (A, /i) is unitarily equivalent to the Hamiltonian obtained through replacing 
g in (3.1) by the dilated function g^ip) :— iJi~^ g{iJi~^p). Introducing the parameter /x is 
thus not so arbitrary as it may seem since ji is in some way related to the width of the 
coupling function. But note that the position of the peak of g is also moved. This is a 
departure from a study which would start with monochromatic photons and enlarge the 
width of their spectrum little by little, without changing the position of the peak. 
This Hamiltonian has invariant subspaces. To describe them and, in the same way, prepare 
the notations for Section IV, let us set 
^!-ad' space spanned by g 

J^n^ := (Ti-ll^^)^ , the space of n photon states, but each photon being in state g. 
Sq := I 0) J^o: and for n > 1, £n '■= | 0) (8) .Fn ® I 1) ® -^n-i) a space we call the 
" n-excitation space" . 

Si^^ := I 0) ® © I 1) J^i'Ji, n > 1, S^^^ = So 



Lemma £n is invariant by H{X, jj), for alln >0 and all ji > 0; En^ C En is invariant 
by H{X,0), for all n > 0. 

As we said when we introduced Notations 1 and 2 in Section II, the method for 
studying the resonances of i/(A,//) is a numerical method in which /j, takes greater and 
greater values, starting from 0. So we begin with giving some properties of H{X, 0). 

111.2 Construction of resonances of H{X, fi) from eigenvalues of H{X, 0). Set- 
ting up 

Proposition 3.1 H{X^ 0) \g(i), the restriction of H{X, 0) to E^^\ has a double infinity 
of eigenvalues 

Co,n(A) = 2-1(1 -Vl + 4nA2) , Ci,n(A) = 2-^(1 + Vl + 4nA2) (3.3) 
For each n, an eigenvector of H{X,0) associated with Ci,n(A) is 

■= {i + nX\^^l{X))-\ I l,^®(-i)) + X Co:^(A) I 0,gn) e 4^) (3.4) 
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Proof £n , two dimensional, is invariant by H{X,0). (3.3) is thus obtained through 
diagonahzing a two-by-two matrix. | 

The choice of the first index in the notation of the eigenvalues is in accordance with 
the principle stated at the end of Section II. By the choice of the second index, we indicate 
that the eigenvector belongs to Sn^ . 

We note that only some part of the double degeneracy mentioned in Section II is 
removed here. Indeed, each Ci,n(A) is still degenerated, since adding an arbitrary number 
of photons orthogonal to g to some state does not change the energy of this state. This 
follows from [ H{X,0), 1 ® c*{h)] — 0, when {g, h) — 0. Let us develop this in order to 
introduce some notations. 

Let ^1,^25- ■• be a basis of functions orthogonal to g. For p > 1, let ^o,p be the 
subspace of £ spanned by 

oo 

0&""-^:=n(l®^*(^^))"^C (3.5) 

i=l 

when n varies from p to infinity; are k non-negative integers, k being arbitrary and 

k 

^ni = n-p. 

In the same way, let Qi^p be the subspace spanned by 

oo 

^(n^,n.,...) -Q {i^c*{g,)r cf>^^l eSn (3.6) 
1=1 

Proposition 3.2 For generic X-values, Qi^p is the eigenspace of H{X, 0) associated with 
the eigenvalue Q^p, for z = 0, 1 and p> 1. 

See hints for the proof in Appendix Bl, which is devoted to the proof of an analogous 
property. 

Let us thus note that eigenstates of H{X, 0) associated with not only are not 
n-photon states but are neither necessarily n-excitation states. This complicates the pic- 
ture of these eigenvalues and of the resonances Zi^n{X, fJ.) they give rise to. 
To be complete, let us still mention that is an eigenvalue of H{X, 0), with, if = n, 
n arbitrary, the associated eigenvectors 

oo 

^^'"''■■■^:=n (3.7) 

1=1 

where f2 denotes the vacuum state of the field. 

Let us now turn to the |U-variation, in order to define and construct resonances 
zo,i(A,/u), 2o,2(A, /u), .. and 2i^i(A, //), zi,2(A,//), .., from their germs ^ij( A) (for the nota- 
tion, see Notation 2 in Section II). 

This construction, as we said, is purely numerical for the moment. It goes step-by-step, 
starting from = 0. We do not look for any existence theorem nor for a complete descrip- 
tion, difficult to get because of the complicated structure of the set of resonances. We just 
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want to obtain numerical values for resonances which have not been considered up to now, 

and might be important. 

This requires a particular choice for g. 

Before that, let us nevertheless give some general indications. 

III. 3 General remarks about the poles of the resolvent of (A, jj) 

In front of so many resonances, it is natural to ask oneself the question: which one is 
the one we are used to, that is to say is there any which could be "associated" with the 
excited state | 1), for the Hamiltonian H{\, 1)? Without going into details here, let us 
explain shortly why the issue is not simple. In particular, given a value Aq of the coupling 
constant, even a small one, why zi^i{\q, /x) is not necessarily the resonance we are used to. 
The resonance we are used to is obtained through restricting H{\, ii) to £i and following 
the resonance which sits at 1 for ^ — 1 and A = 0. We follow it as A increases from to Aq. 
Since its position at A = does not depend on /i, it amounts to following the resonance 
along the path of Figure la, from its value 1 at the origin of the path. 



1 [ 

X.Q X. X. 

(a) (b) 
Figure 1 . Two paths enabling to construct two non-necessarily identical resonances. 

It is the usual perturbative approach, in which there is no need to introduce /U. On the 
contrary, -2i,i(Ao, 1), also defined with £i (see III. 4. 1.2), is the limit of -2i,i(Ao, /u), when 
increases from (then -2i,i(Ao, 0) = Ci.il-^o)) to 1. Since Ci,i(Ao), by definition, tends to 1 
when Ao tends to 0, -2i,i(Ao, 1) is the resonance obtained through following the resonance 
along the path of Figure lb, from the value 1 at the origin of this path. Now, in some 
examples, functions 2o.i(A, /u) and -2i.i(A, /u) are two branches of a unique analytic 2-variable 
function, which, even when restricted to IR^, has branch points. These real branch points 
may lie inside the rectangle [0, Aq] x [0, 1] of the (A, /x) plane, so that the two previous values 
of the resonances at the common end (Aq, 1) of the two previous paths may be different. 
The A and jj, variations do not necessarily commute"'^^. (Let us note here that it is thus 
dangerous to move A and ^ simultaneously without caution, as we did it for instance in 
the second paragraph of Section III of Ref. 17). This phenomenon is at the origin of the 
remark in parentheses in the first paragraph of Section III. 

This seems to be a drawback of the introduction of parameter ji. But this introduction 
enables us to see resonances for the physical Hamiltonian {fi = 1) which the perturbative 
approach does not give so easily. Indeed, when A is the only parameter and when it is 
varied from Aq to 0, the continuous variation of the resonance siting at -2i,i(Ao, /u) for A = Aq 
may lead to 1, but also to infinity, or to a pole of g (see Figure 2, where the pole is —i). 
The result depends on ji and g. But these two latter limits are too singular points to start 
a calculation in their neighbourhood. 
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Since the difficulties of a general study are due to branch points whose positions 
depend on we choose a particular g. We now come back to the study of the announced 
particular case. 

III. 4 Construction of resonances of i?(A,//) from eigenvalues of H{X, 0). Cal- 
culations for a particular g. 

We fix the coupling constant at 0.1 and take the particular function we already used-'^°: 

This is a simple rational function which exhibits the type of singularity that actual coupling 
functions may have. (See for instance matrix elements of the interaction Hamiltonian for 
hydrogenic atoms in the electromagnetic field^^) 

In Section III. 4.1, we recall the definition of 2o,i (-^5 A*) a-iid i (A, /x) and the known formulas 
through which they are obtained. We give their values for various ^ and A. We use them 
later on. Then, in Section III. 4. 2, we define resonances 2o.2(A, /u) and 2i^2(A, //) and give 
some approximate values. An evaluation of the errors is given in Appendix A. 

III.4.1 Brief review about two resonances obtained with only one photon: 

2:0,1(0.1, /i) and 2:1,1(0.1, /x) 

They are poles of (1 f2 | [z — H \£^]~'^ \ 1 (8) f2) or of its analytic continuation. 

III. 4. 1.1 Resonance ^o,i(0.1,/u) 

2:^,1(0.1, h) are zeros of /(A, .), defined by 

f(X,H,z):^z-l-2X^ f -^^^dp (3.9) 

Jo z- lip 

in the cut plane C/IR"*", or of its analytic continuation through the cut, clockwise^ 
For ji < //c(A) := 2A^ ^^^dp, /(A, .) has only one zero; it is on the negative real axis 
and is if and only if ^ = /^c(A). When ji tends to 0, it tends to Co, 1 (0.1) and therefore 
we denote it by 2:0,1 (0.1, /u). We recall that the corresponding normed eigenvector, which 

we denote by 00*^1 1'^' 1A1 proportional to 

<J(A,^):= |l,O) + A|0, — -) (3.10) 

2o,i(A,/i) -/ibi 

For > //c(A), -20,1 (A, .) is defined as a zero of the analytic continuation 

/+(A, 2) 2 - 1 - 2A2 / ^^dp + ^m-g{-f (3.11) 

Jo z- lip II II 

of f{X,iJ,,.), clockwise through the cut. When A*c(A), it connects to the values for 
li < Hc{X). For H varying from to 2, some values of 2:0,1 (0.1, h) are given in Table 1 and 
Figure 1. 



10^ li 





0.1 


1 


3 


6 


6.2 


6.36 


6.366 


10* 20.1(0.1, ;U) 


-99 


-94 


-68 


-34 


-2.5 


-1.1 


0.04 


0.001 



Table 1. Values of 2o,i(0.1, /x), for some values of 11 below /ic(O.l). 
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Figure 2. The resonance zo,i(0.1, for n e [0,2], first in H , then in the 
second sheet of the complex plane 

The physical value for ^ = 1 is 0.11 - 0.95 i. For /i = 0, 2o,i(0-l'0) = Co,i(0-l) ~ 
-0.0099 ; the graph goes through for ~ 6.36 lO'^ and 2o,i(0-1. 2) = 0.13 - 1.97 i. 

In order to connect these results to the usual perturbative treatment, we drew a 
dashed line in Figure 2. It describes the continuous move of the resonance which sits at 
the point 2:0,1 (0.1, 1) for A = 0.1, when A decreases from this value to 0. One can see that 
the limit is not 1, the energy of the naked excited state, but —i, a pole of g. This makes a 
difference with what occurs for 2:1,1 (A, 1), as we will see it just below. To distinguish the 
two behaviours, we called a resonance which, as a function of the only A variable, does 
not tend to the excited state energy when A tends to a "non standard" resonance-^^. Be 
careful that a resonance may be standard for some /i's and non standard for others. 

III.4.1.2 Resonance 21 i(0.1,/i) 

It is another zero of /+(A, /U, .), which tends to Ci,i(0.1), when — > 0. Its displacement 
when n varies is given by the full line of Figure 3. 

XM = 






995 / 
/ 


1 . 005 


1 . 0^ 


= 2 


/ 

/ / 
/ / 

1 / 

/ 

\ / 

\ 76. 008 
\ / 






n = 


1 ^^-0.01 
-0^0^ 







Figure 3. zi^i{0.1, fj,), in the second sheet, for fj, e [0, 2] (full line) 

We have 2i,i(0.1,0) = 1.0099 , 2i,i(0.1,l) = 0.997 - 0.010 i , 2i,i(0.1,2) = 0.995 - 
0.0032 i. It can be shown that the resonance which sits at point 21,1 (A, 1) for A = 0.1 
moves continuously to 1 along the dashed line of Figure 3 when A goes to 0. Therefore, 
referring to the discussion of Section III. 3, we can say that 2:i,i(0.1, 1) is the resonance 
usually associated with the excited state. 
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We now come to the original part of this Section III. 

III.4.2 Two other resonances, obtained with two photons 

Among the various resonances tending to Ci,2(A) when ^ goes to 0, we select poles 
of matrix elements of [z — H{X, fi) ffj]"^' '^^^ their analytic continuation. Note 

that another choice among those mentioned after Notation 2 would have been to look for 
poles of matrix elements of [z — H{\,ii) ffg]""^, since we have eigenvectors of H{X,Qi) in 
£3 associated with the same eigenvalue Ci,2(A): for instance (l ® c* {g\))(j)i^2- But we make 
no claim to being complete. We just want to give an example of two resonances which are 
not considered usually. 

Two new difficulties now appear. The first one is that, contrary to the £i-sector case, 
resonances in the £^2-sector are given by zeros of a function D(A,//, .) which is no longer 
explicit, since we will see that it is the sum of a series. Up to now, we can only calculate 
these zeros approximately, through cutting the series after the first non trivial term. The 
zeros of this truncated function are the approximate values we consider, for the resonances 
in this sector. 

The second difficulty is that this truncated function has several branch points; as a con- 
sequence, it has several analytic continuations. There is no reason why it should not be 
the same for D{\,fi, .) itself, although it is possible. Therefore, when considering a zero, 
one must tell which branch is in question; we give a notation for the zeros of the various 
branches later on. 

III.4.2.1 These resonances as zeros of a function D{\ji, .) 

In £21 looking for the eigenvalues leads to the following proposition. 

Proposition 3.3 (Eigenvalues of H{\, fi) as zeros of a multivalued function) A 
being fixed and n and z being two parameters satisfying /i > and z e IR~ , let V^^z be the 
Predholm function (see Ref. 22, p. 68j associated with the integral equation 



with 



(p) j K^,,{p, q) ifx,M dq = (3.12) 



where 



,AP,q) (^_i_^|^|_A2 (T^,,<7)(p))(^-MN + |g|)) ^ ^ 
i-.M^^-IJ^,'^' (3.") 

Let z — > D{X, jd, z) be the multivalued function which, for z e 1R~, equals ©^^^(A^). 

In the case where one of the zeros of D{X, ^, .) is real negative, this zero is an eigenvalue 
of H{X, jj). Let us denote it by ^o,2(A, jj). The associated eigenvector is in £2, proportional 



to 



<^(A..):^ |l..U + v^A|0.^-^^^^|^^^-^) (3.15) 



where <^a,/u « solution of 

<P\Ap) J ^m,$o,2(a,m)(p>?) 'PxAq) dq = (3.16) 
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Proof I l,{px^^) + I 0, XA,^t) is an eigenvector of H{X,ij,) associated with the 

eigenvalue z if and only if (3.12) holds and X\,fj, = X{z — /u(|pi| + \p2\)) 9 V V'a,/x- 
According to Fredholm's theory, (3.12) has a non trivial solution only if ^^(A^) = 0. 
The Proposition follows from the fact that (3.14) and (3.13) are defined if 2; is a negative 
number. | 

Remark If /x is set to in (3.12), (3.13) et (3.14), one finds that (3.12) implies z = Ci,2(A). 
If the continuity with respect to /x could be proved, we would get that the limit of ^0,2 (A, //) 
when fj, goes to is Co, 2 (A). We do not know how to prove this at the moment, since we do 
not know all the zeros of -D(A, /U, .). Nevertheless, the result of the approximate calculation 
of Section III. 4. 2. 4 is along this line. That is why we eventually change notation ^o,2(A,//) 
for Zq 2 (A, /x), according to Notation 2. (The upper index is explained later on.) 

To switch from eigenvalues to resonances, let us take h G L^(IR^), z s.t. $52; > and 
ip := I 1,/i). Let us introduce H^^ziPiQ.i^'^)-, the resolvent kernel of equation (3.12) (see 
Ref. 22, p. 63). It can be shown that 

(V^, [z - H{X, /.)]-! V) = 1 - A' y" H^Ap, <1, >^^)h{p) Kq) dpdq (3.17) 

We know (Ref. 22, p. 58 and 63) that the only singular points of H^^p^q^ .) are the 
solutions of V^^X) = 0. The zeros of analytic continuations of function L>(A,//,.) of 
Proposition 3.3 will thus give us poles of the left-hand side of (3.17), that is to say reso- 
nances. 

The calculation will be an approximate one. The result for = 1 is given by lines two and 
four of Table 3, in Section III. 5. 

III.4.2.2 Zeros of D(A, .) approached by zeros of a function D^^\X, /i, .) 

~" i-iy 

nl 



Proposition 3.4 For z < 0, -D(A, /i, z) = ^ j — C'n(A, z) with Cq = 1, 

n=0 

Cn{X, fi, z) = X'^'^ J Dn{lJ^,Z,pi,- ■ ■ ,Pn) dpi- ■ -dpn (3.18) 

and, for 1 < n < 3, 



i=n 



Proof The D^s are determinants given by Fredholm's theory. We calculated them for 
n < 3; the result is given in (3.19). | 

The analytic structure (poles, branch points) of -D(A,/i,.) is difficult to determine. 
Indeed, £)(A, /i, .) is a series of terms each of which has a different analytic structure. 
Moreover, this structure is not simple. Therefore we will limit ourselves with replacing the 
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search for zeros of D{X,fi, z) by the search for the zeros of the sum of the first two terms 
of the series: 

D^^\X, fi, z) := 1 - Ci{X, 1^, z) (3.20) 

The first corrections to this approximation are discussed in Appendix A. From Proposition 
3.3, we get, if 7^ 0, 



where 



i?W(A,/x,^) = l 



V'(A,//,^,q) := 



2A^ 



1^ Jo 



ip{X,l^,z,q) dq 



(3.21) 



(3.21') 



{z-2\q\) f{X,^i,z-\q\) 

andD(^)iXAz)= %-X'S - 

III.4.2.3 Analytic properties of D^^^ 

Proposition 3.5 For 7^ 0, D^^\X, .) has at least three branch points: 0, 2:0,1 (A, fi) 
and 2i,i(A, n). 

Proof Let us recall that /(A, .) is analytic in the complex plane cut along IR^, with 
a branch point at 0. The analytic continuation from the upper half- plane has a pole at 
—ifi. An explicit expression of / in the upper half-plane is, for z 7^ ifj,, 



2 2^^+^^7TZ + 2^Z 

f{X,ii,z) = z-l + X 



^ - 7iz^ + 4/12;^ ( log /X — log(— 2)) 



(3.22) 



The continuation clockwise across the cut, which we denote by /_|_(A, //, 2), is obtained 
through adding Am^g{ii~^z)'^ to the above expression. It is convenient to rather introduce 

the function f{X,iJ,,z) which coincides with f{X,fj,,z) in the upper half-plane and has a 
cut along ilR~ . 

For ^z > 0, the poles of ip{X, yu, z, .) (see (3.21')) in the complex plane cut along z + iJR'^ 
are ±z//, z/2 and qi{z, fx, A) := z — Zi^i (A, fx). Depending on whether fx is smaller or greater 
than //c(A), the 2;-dependent poles sit at places schematically shown in Figure 4. 



qi(z,n) 



z/2 



z/2 



4a 



4b 



Figure 4. Branch point and poles of \^(k,\i, z,. ) , defined by (3.21') 



12 



This follows from the position of Zi^i{X,iJ,), given by the curves in Figures 2 and 3. 

If /i = /ic(A), the pole qo{z, coincides with z. (In Figure 4, z has a positive real part, 

but it could be negative, as well.) 

When z enters the lower half-plane along a path 7 (for instance the dotted lines in Figures 
7 and 8), the integration path in (3.21) may have to be deformed in order to be kept away 
from some of the three poles or from the branch point z of il}{X.i 2;, .). Two different paths 
7 and 7' will not necessarily yield the same result. For instance, if 7 crosses Zi^i (A, fi) +IR''", 
qi{Zj n) crosses IR"*"; the integration path has thus to be deformed, whereas it is not the 
case if 7 crosses ;2i,i(A, /x) — IR'*'. As a consequence, 2;o,i(A,/x) and ;2i,i(A,/x) are branch 
points. 

If z comes to zo^i{X, fx) (resp. zi^i{X,ii)), then the pole qo{z,ii) (resp. qi{z,ii)) comes to 
and the integral in (3.21) is singular. 

To describe the various branches of D^^^ readily, we need a notation for the homotopy 
classes of paths in Xx^^ :— (C\ {0, zq^i{X, /x), zi^i{X, /x)}. (Two paths are homotopic if they 
can be continuously deformed into one another, in -^a,/u-) Depending on whether /x is 
smaller or greater than fic{X), the relative position of the three branch points is different. 
We refer to Figure 5 for the definition of the fundamental class of paths a, ao, ai. The base 
point B is chosen real and smaller than Co, 2 (A). 




Figure 5. Fundamental paths in the z-plane. 

The path 'jo{z) is defined as the polyline going through points B + ze, ^{z) + ie and 
ending at z, for an arbitrary e > 0. Every path 7 from B to z can be expressed by means 
of O, Qq, tti and 7o(-2)- It defines a homotopy class [70(2;) ^7] G 7Vi{Xx^^). Conversely, with 
each class [/] of tvi{Xx,ij,), we can associate the homotopy class of paths ^o{z)l going from 
B to z. With each pair consisting of [I] in 7ri{Xx,ij,) and z, we can associate the analytic 
continuation in z of D^^^/i, .) along the path ^o{z)l. We denote the value at the end of 
the path by D'^^^{fj,, z). From now on, we will omit to mention the A variable (which has 

been fixed to 0.1.) We denote by z^''\fj,) a point such that D^^^ (^fj,, z^''\ii)) — 0. When /i 

varies, we assume that this point varies continuously. It is denoted by zf\{ii) if its limit 
when fx goes to is Ci- 

These zf\{ii) are approximate values for the resonances we are considering. As regards 
their physical meaning, we refer to the short comment in Section III. 5. 

We now give some values of these functions for various values. 
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III.4.2.4 Values of resonance -Zq 2(a*) 
(i) II < 6.3662 10-3 

(a) [I] — 1. For n close to 0, we look for a zero of D^^\p,.) in the neifflibourliood 
of Co,2 = -1-96 10-2. .) is well defined by (3.21) in the neighbourhood of every 

negative real number. A calculation on a computer yields negative real zeros of this 
expression, for ji small (see Table 2, first column). When ji increases up to a certain value 
Ijl'^21 close to lie '■= A*c(10-^), the same formula still gives a negative real zero. The first 
column of Table 2 gives its values for < < 6.3662 10-^ ~ jic- 



lOV 















-196 











103.9 


0.1 


-185.5 


1.65 10-3 


3.95 10-6 


1.04 10-^ 


105.3 


1 


-135.7 


13.14 10-3 


10-^ 


6.26 10-^ 


116 


3 


-69.8 


0.038 


7.93 10-^ 


10-5 


143.8 


6 


-9.3 


0.099 


7.7 10-3 


3.68 10-4 


247.1 


6.2 


-6.6 


0.106 


9.3 10-3 


5.2 lO-'^ 


265.3 


6.36 


-4.4 


0.111 


1.09 10-2 


7.1 10-^ 


284.2 


6.3662 


-4.328 


0.112 


1.10 10-2 


7.17 10-4 


285 



Table 2. Values of -20,2 ( A*); some values of ji below jic- 

Through comparing Table 2 with Table 1, one sees that this zero is smaller than ^0,1 (a*)- 
In accordance with the notations of the end of Section III. 4. 2. 3, and with the upper index 
1 denoting the unity element in 7ri(X^), this zero may then be denoted by Zq 2(1^) ^ since 
the branch point 2o,i(a*) does not belong to the interval Bzq2{iA- Figure 6a shows the 
three cuts of D^^\ii, .) for 11 < 6.3662 10-3 and this zero Zq 2{iA- 

It has not been possible to determine the place of //^2 ^i^^ respect to /ic(10-"'^). We 
will come back to this point in Appendix A when we estimate the errors made in the 
approximation. 

(h) [/] = ao-' a zero associated with another branch. Figure 6a also schematically 
shows another zero: z^q'^[ii). The notation indicates that this value is now the zero of the 

continuation of D'^^^p, .) along the path 7' := 'Jq(^Zq^2\p)) '^Oi ^ path which goes around 
the branch point zo^i{fj,). Figure 6b shows the two paths followed by qo{z, //) when z goes 

along 70 (4^2 (/^)) 7'. 



z (n) 

0,2 



zo,l(^l) 




z - plane 



z 1 1 (1^) 



Fig. 6a 



q - plane 



q( z' (H) ,n) 

o 0,2 




> .. . • 



q( z^ (n),n) Fig. 6b 

0,2 



[11 -3 

Figure 6. Paths defining resonances z^ ^ for |x < 6 10 
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This 2;o°2(/i), also close to Co, 2 for fJ> ^ 0, is a zero of D^J^/j,, .) which, for 3^(2;) < (and 
also for' 3^(2;) > 0,$J(2;) > 0), reads 



{z-2\q\) f{ix,z-\q\) 



dq + 4z7r 



g{qo{z,fi)/fiy 



H {z -2qo{z,ii))dJ{zo^i{iJ,X),n,X) 



(3.23) 



with A = 0.1. This zero is no longer real. 

(ii) For |U ~ things are not clear. 

(Hi) For // > 7 10~^, the branch point 20,1 of D^^^ is in {z; ^(z) > 0, ^{z) < 0} (see 
Figure 2); it is shown in Figure 7. Among the various analytic continuations of D^^^ across 
H"*", we consider d[^\ Expressions are given below (formula (3.24) and the two last lines 
of the section). The zero of d[^\ denoted by -20,2 (a*)) follows the curve of Figure 7 when /j, 
varies from 7 10~^ to 1. For = 1, the way D^^^ is analytically continued to the zero of 
d\^^ is shown by the dotted line. 



^/i = 0.007 

ro5 



. 15 




Figure 7. -20,2 (a*)) i^i the complex plane, for /j, e [0.007, 1] 

For = 7 10-3, 4 2()") = 2.8 10"^ - 2.4 10"^ i, close to 0. Values of ^0,2 (/^) for between 
6.37 10~3 and 7 10"^ are difficult to get. Actually, their determination is not useful 
since Zq 2(1^) is only an approximate value of the exact resonance. We have 42(0-1,1) = 
0.216-1.9 i. For ^z < 0, ^{z) > and > 7.35 lO'^, stiU with A = 0.1, 



fJ' ^ J[o,U{z)]u[U{z),z] - 20-) /+(//, z-q) 

^[0,SR(«)]u[SR(z),^] (^ - ^Q) JKfJ'i^- Q) ^fJ' 

2iV(i^fMM) (2zo,i(/^)-^)-' dj+ii,, zoM)-' 
A* 



(3.24) 



This formula is derived from (3.21) in the following way: when z enters the quadrant 
{2;; 3^(2;) > 0,3^(2;) < 0}, the cut in Figure 4b drags the integration contour along, which 
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yields the first two terms. The last two terms come from the residues of the poles 2;/2 and 
qo{z, IJ,) which cross IR"*" {qi{z, jj,) does not cross M''"). For 10^// G [7,7.3], the expression 
of D^\iJ,j z) does not contain the residue term at go since 7o(^o2(a*)) does not cross 
zo,M+M+. 

III. 4. 2. 5 Values of resonance 2;} 2(a*) 

It can be shown numerically that the zero of d\^\h, .) which tends to ^1^2 when n 
tends to 0, denoted by 2i^2(a*); is a zero in {z; ?R.{z) > 0, ^{z) < 0} of A(/i, .), where 

A(/x,^):=l- 



fJ' ^ i[o,SR(«)]u[SR(z),2] - 2g) f+ifJ', z- q) 



2ing'{^^^hM) {2z,,M - z)-' ^M^^,z,M)-') (3.25) 



^{n^z) is the expression of D^^\ijl,z) in the neighbourhood of the considered zero, but 
not everywhere in the lower half-plane; for example, these two functions differ at a point 
z such that '^{z) < 9(2;o,i(//)). The variation of -2i^2(a*) with /j, is given by Figure 8. 





11 = 




-0.5 


1.03 ~TTCW — 


— i-,-ai_^.0 6 
i 


-1 






-1.5 






-2 


^11 = 2 





Figure 8. zX ^ip)-, in the complex plane, for e [3 10 ^, 2] 



This resonance starts from 2~^(1 + \J\ + 8A^) = 1.01962 for = and goes through 
1.043 — 1.127 z for = 1. We note that this value is much farther from the real axis than 
the resonance 2;i^i(0.1, 1) of sector S\. As in Figure 7, the dotted line in Figure 8 indicates 

the path of the analytic continuation of D'^^^ to the zero (of d!^^), for // = 1. The origm 
is a branch point. The other two branch points, 2o,i(0.1, 1) and Zi^i(0.1, 1), whose values 
are recalled in Table 3, cannot be drawn at the scale on the figure. 
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III.5 Conclusion of Section III 

Table 3 gathers the four results we obtained for A = 0.1 and = 1. It gives four 
resonances among all those of the physical (i.e = 1) Hamiltonian, when the function g 
is given by (3.8). 




0.13- 1.97 z 
0.216-1.9 i 
0.997- 0.010 i 
1.043 - 1.127 i 



Table 3. Four resonances (approximate values for two of them), 

for A = 0.1 and /j, — 1 

We see that the usual resonance zi^i is the only one near the real axis. However, for small 
values of this is no longer true, as it can be seen from Figures 2, 7 and 8, or Table 
2. Therefore, resonances for small fi may play an important role. Now, in the harmonic- 
oscillator example mentioned in the introduction, what plays the role of parameter fj. is 5~^, 
the inverse of the spatial extension of the states. Thus, for extended states, /x is small, and 
we recover our motivation for the study of all the resonances of atom- field Hamiltonians. 

One should of course discuss which of these resonances have a physical meaning and 
make this meaning precise. One should also find actual situations in which these resonances 
or eigenvalues can be seen readily. Unfortunately, some more work has to be done. We 
think that the physical meaning should appear quite easily in the case of eigenvectors of 
the Hamiltonian. In £i (see (3.10), Table 1 and Figure 2), we have an example of such a 
state. As regards the restriction of H{\ji) to £21 we did found approximate real values 
for the resonances. But we have not shown that (real) eigenvalues do exist. The study has 
to be carried on. Now, as pure resonances are concerned, we think that their study cannot 
simply be an academic question, since eigenvectors change continuously into resonances 
when parameters of the physical system are varied. 

In any case, the 2-level model is perhaps too simple to find a concrete application. 
The next section is a step towards more realistic models. 

IV Hybrid states for a three-level atom coupled to photons 

To be in a position to describe mixed states in more realistic models, we not only must 
be able to consider several photons, as in Section III, but we must also be able to consider 
several atomic or molecular levels. The present section is a preliminary study devoted to 
mixed states when S is three-level atom coupled to the radiation by a Hamiltonian of 
type (2.2). 

The reason why the study is only a preliminary one is that there come some additional 
difficulty, together with those already mentioned in Section III: for /i > 0, none of the spaces 
with a bounded excitation number is stable by the evolution operator, contrary to the two- 
level case. This is due to the possibility of the transition | 0) -|- photon — > | 2) ( | 2) is the 
second excited state), through which the excitation number may increase. For this reason, 
the determination of the Hamiltonian eigenvalues is already not simple for fi = 0. Since we 
assume that the structure of the resonance set will roughly be conserved when fi takes non 
zero values, we must study the (/x = 0)-problem first. This is the subject of this Section IV. 
The displacement of these resonances when /j, becomes non-zero will not be examined in 
the paper. 
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A typical result is illustrated by Figure 12. 
IV. 1 Notations and Hamiltonian 

The atom has three levels with energies cq = 0, ei and 62, corresponding to states 
I 0), I 1) and I 2). /oi,/i2,/o2 being three normed functions in L^(IR), the Hamiltonian 
is 

if(A,/i) :=ifo(A,/x) + Ao2 V (4.1) 

with 

i?o(A, y) := (ei I 1)(1 I +62 I 2)(2 I ) ® 1 + ;u 1 ® i?rad+ 

Aoi( I 1)(0 I ® c(/oi)+ I 0)(1 I ® c*(/oi)) +Ai2( I 2)(1 I ® c(/i2)+ | 1)(2 | ® c*(/i2)) (4.2) 

V :=| 2)(0 I ® c(/o2)+ I 0)(2 | ® c*(/o2) (4.3) 

-f^rad is as in Section III.l. Transition | 0) I 2) is distinguished from | 0) ^ | 1) and 
I 1) — > I 2) for a physical reason and also to prepare a perturbative calculus. 

We use the following notations, modeled on those of Section III.l. As regards the 
photons, we denote the Fock space by the n-photon space by !Fni the vacuum state by 
ri, and we set 

^rad- space spanned by /oi, /12, J^n^ '■= ('^rad)^ • photon-state space, but 

each photon being restricted to be in "^^^ad -^^^^ '■~® ■ 

n 

^rad" space Spanned by /oi, /12, /025 '■= ("^rad) ( ^^^^ photon is restricted 
to be in ni%) and .F^^) :=® . 

n 

As regards the atom-photon system, Sn and £ have been defined in Section III.l. We set 

£:f ) := I 0) ® JP-{^^ ® k = 2,3 

£^n^ := I 0) ® JSi''^ e I 1) (8) r^nl^ ® I 2) (g) JF^^Ja, = 2, 3, n > 2 

£:(^) := ® 

n=0 

A := (Aoi, A12, A02) 
We also introduce sq (/oi,/o2), si := (/oi,/i2), ^2 := (/o2,/i2). 

We will use the letter (f) to indicate eigenvectors of ifo(A,0); a priori, they depend on Aqi 
and A12. We will use x to indicate eigenvectors of -ff(A, 0); they also depend on Ao2- 

We aim at getting the eigenvalues of i^(A, 0), an operator which we simply write H{X), 
from now on. When the variable /x is not mentioned, it will be assumed to be 0. 
To this end, we take up the idea mentioned in Section II consisting in perturbing H{0) 
through introducing the interaction step by step. The doubly infinite degeneracy of the 
eigenvalues cq, ei and 62 of H{0), due to the arbitrariness of the number of photons and 
the arbitrariness of the state of each photon in the corresponding eigenvectors, is partially 
removed at each step. The first perturbation will be the addition of the Aqi and A12 terms 
of (4.2) to H{0). It is described in Section IV. 2. The second perturbation will be the 
supplementary addition of A02 V. It is described in Section IV. 3. 
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IV. 2 Perturbation with respect to Aqi and A12. First splitting of cq, ei and 



62 



Here we are interested in Hq{X) := Ho{X,0); the interaction X02V is switched off. 
First, in Sq, \ 0,0), is an eigenvector associated with the eigenvalue 0. 

(2) 

IV. 2.1 Three eigenvectors of Ho{X) in the 1-excitation space and the 

three associated eigenvalues 

(2) 

The space , three-dimensional, is invariant by Ho{X). It is the direct sum of the 
eigen-subspaces | 1) <S> Tq'^ and | 0) (g) T^'^ of -ffo(O), associated with the eigenvalues ei 
and Co respectively, and dim. = 1 and dim. J^[^^ = 2. The first perturbation will shift 
ei and split cq into two eigenvalues, as it is represented in the first two columns of Figure 
9. 

62 




C: 



1,1 



1 



(Co,l)2 

(Co,i)i 



l{0) 
'1,1 



(<J)2 



|1,0>, |0,1> 



|0,1> 



(<^0,l)l |1,0), |0,1) 



levels of the atom 
decoupled from the field 



1-excitation eigen- 
mixed levels values 



dimension of the field 
state-space (apart eigenvectors and 

from degeneracy) their components ^""^^ 

(2^ 

Figure 9. Levels associated with eigenvectors of Hq{X) in , 
therefore without any spectator-photon 



The third column gives the notation of the perturbed eigenvalues. In the fourth 
column, we recalled the above mentioned number 1 and 2, which are the dimensions of 
the projections of Sl^"^ on | z) (S> JF^^). They also are the sum of the dimensions of the 
eigenspaces of H{X) L(2) associated with eigenvalues which tend to when A tends to 0. 

(2) 

"Apart from degeneracy" indicates that in considering we do not take account of the 
additional degeneracy of the eigenvalues due to photons in states orthogonal to /oi and 
/12; we call these photons spectator- photons. This degeneracy is explained in Proposition 
4.2 and is removed by the second perturbation studied in Section IV. 3. 
The eigenvalues obtained through the first perturbation and the associated eigenvectors 

are given by the following proposition. 

(2) 

Proposition 4.1 In Si , i//oi 7^ fi2, Ho{X) has three eigenvalues: 

(Co,i)i(Aoi, A12) = 2-i(ei - (4-4) 

(Co,i)2(^oi,Ai2) = (4.5) 

\ i, n) schematically indicates a 7l-photon state with the atom in state i . 
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Ci,i(V, Ai2) = 2-1(61 + (4.6) 
The associated normed eigenvectors are 

(<J)i(Aoi, A12) = (1 + (Co,i)r'A^i)"'/'( \l,n) + (Co,i)r'Aoi I 0,/oi)) (4.7) 

«J)2(Aoi, A12) = (1 - \si\')-'/' I 0, A2 - si/01). (4.8) 

<;(Aoi, A12) = (1 + (Ci,i)-'A^i)"'/'( I 1,0) + (Ci,i)-'Aoi I 0,/oi)) (4.9) 
The calculus is straightforward since the dimension of is three. | 



To the second order with respect to A, we get (Co, 1)1 = —-^Aoi and = ei + ^Aqi. 
Note that the perturbed eigenvalues (4.4) to (4.6) do not depend on A12. They tend to 
0, and ei, respectively, if Aqi tends to 0. The notation is thus in accordance with 
the rules stated in Section II, as far as the first index is concerned; index i refers to the 
unperturbed level in the following way: limAoi^o(Co,i)i — , limAm^oCi,! = ^i- To 
distinguish perturbed eigenvalues which tend to the same if (Aoi,Ai2) (0,0), we 
chose to put the index p, indicating that the eigenvectors are in Sp (as in Section III. 2). 
Since several eigenvalues with the same indices i and p may still have the same limit when 
(Aoi, A12) tends to (0, 0), an additional index j is used to number them. 

We have (^{,°J) ^(Aqi, A12) - \ 0, /oi) and (Aqi, A12) | 1,0), if Aqi tends to 0. 

If foi = /i2 7 the dimension of s[^^ is two and the eigenvector {4>^\)2 disappears. 

As we said before, these eigenvalues are actually infinitely degenerated; indeed, adding 
photons orthogonal to {/oi, /12} gives an eigenstate with the same energy, since the energy 
of the photons is not taken into account in Ho{X). Let us state this fact precisely, with 
some notations which will be useful later on. 

IV. 2. 2 Other eigenvectors of Ho{X) in the (n+l)-excitation space Sn+i, n > 0, 
associated with the same eigenvalues. 

Proposition 4.2 Let gi, g2, - ■ ■ be an orthonormal basis of functions orthogonal to /oi 
and /12. 

(i) Let ^0,1,1 be the subspace of £ spanned by the normed vectors 

00 

i<f>t+r-\ (^^0-Ml®c*(^,))"^ (4.10) 

where the rii 's are k non-negative integers and k is arbitrary. These vectors are in £n+i 

k 

ifn — '^^ni. For generic values o/Aqi and X12, Go,i,i is the eigenspace of Hq{X) associated 

i=l 

with eigenvalue (Co,i)i- 

(ii) If foi 7^ f 12, let us set go := (1 — |sip)~2 (/12 — si/01) and let /qI be the subspace 
of Ti-r ad orthogonal to foi , spanned by the gi's, i = 0, 1, ■ ■ ■. The eigenspace of Ho{X) 
associated with eigenvalue (Co,i)2 is | 0) ® J^{foi), where J^{foi) is the Fock space built 
with foi- We set 

'•••^)2 (<1)2 (4-11) 

i=0 
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(in) Lastly, let Qi^i he the subspace of £ spanned by the normed vectors 



/ (ni,n2,---) 
h,n+l 



(4.12) 



For generic values o/Aqi and A12, Qi 1 is the eigenspace of Hq{X) associated with eigenvalue 

Ci.i- 

Thatn,=i {l^c*{g,)r (<l)i, n.=o {^^c^ig,)^ (<1)2 and a=i (l®c*(^,))"^ ^S^J 
are eigenvectors follows from the fact that [1 <S> c* {g),Ho{X)] = if is orthogonal to H^^^. 
In Appendix Bl, we explain how the whole eigenspace can be determined, for each of the 
three eigenvalues. | 

Let g (which is no longer function (3.8)) be a linear combination of the gi^s. Propo- 
sition 4.2 may be stated in the following terms: the three eigenvalues are twice infinitely 
degenerated; firstly through the number of spectator-photons (the n variable), and sec- 
ondly through the infinity of possible states for each spectator- photon ( the g variable). In 
the level diagrams, if we symbolise the degeneracy of an eigenvalue (due to the possibility 
of one spectator-photon in the eigenvector) by a dotted line, then in the case where at 
most one spectator-photon is present, the degeneracy of the levels may be represented in 
the following way: 




Ci,i 



(Co. 
(Co^ 



2JV 



2,0) 



two-excitation 
mixed levels, 

levels of the decoupled with spectator- eigen dimension of the 

atom photons^"*^^ values field state-space eigenvectors 

Figure 10. Levels of Figure 9, with eigenvectors of Hq{\) in £^2; 
at most one spectator-photon 
The degeneracy of the eigenvalue is greater if we admit eigenvectors with a greater 
number of spectator-photons, with all their possible states. 

But considering a total number of photons greater than one yields other eigenvalues. 
For instance, wc claim that eigenvectors in 82 are not necessarily of the form given in the 
last column of Figure 10. We are going to see that the set of perturbed eigenvalues coming 
from one given unperturbed energy changes when the number of photons coupled to 
the atom changes. This was already the case for the two-level atom of Section III (see for 

dotted lines are slightly shifted from the level for reading purpose 
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instance Proposition 3.1 of Section III. 2, for /i = 0, where the eigenvalue depended on the 
space in which the eigenvector was looked for). 

IV. 2. 3 Other eigenvectors of Ho{X) in £n+i and other eigenvalues 

In this section, we assume /oi ^ fi2- 

Proposition 4.3 Ho{X) has an infinity of eigenvalues different from (Co,i)i: 
(Co,i)2 et Ci,i. 

Proof £n is invariant by Hq(X). The subspace £n^ of Sn is also invariant. 

(2) 

Let us first consider S2 , six-dimensional. If the coupling constants are small, there are 
six eigenvalues: when Aqi and A12 tend to 0, one of them tends to 62 (it is denoted by 
(2,2) , two tend to ei (denoted by (Ci,2)i and (Ci,2)2) and three, one of which is zero, tend 
to (denoted by (Co, 2)1 = 0, (Co, 2)2 and (Co, 2)3)- They are obtained through diagonalizing 
a six by six matrix (see Appendix B2). The calculation is straightforward although the 
result has not a simple expression. In this 2-excitation space, that two eigenvalues tend to 
ei is due to the fact that there are two possible photon states, and that three eigenvalues 
tend to is due to the fact that there are three possible independent states for the two 
photons. The fourth column of Figure 11 recall these numbers. The eigenvector associated 
with (Co, 2)1 = is I 0,gfo V go); we met it under the form {4>Q^f' ^)2 (sec (4.11)). 
We already found other eigenvectors in the same space 82- They were built from eigenvec- 
tors in £[^\ They are for instance {(l>Q'f' ^)i and \ with notations of Proposition 

4.2; they are in S2 if (71 G ^rad- '^^^ corresponding eigenvalues, (Co, 1)1 and Ci,i, are 
different from the six we just saw (see Appendix B2), except possibly for particular values 
of A. 

Since the notations are a bit heavy, we again represent the levels in Figure 11. 

|2,0), |0,2) 




(Cl,2)2 
(Cl,2)l 



(Co, 2)3 

(Co,2)l 
(Co,2)2 



(<l)3 
(<l)l 



eigenvectors 

and their contents 



levels of the dimension of the 

decoupled 2-excitation eigen- field state-space 

atom mixed levels values (apart from deg.) 

Figure 11. Levels associated with eigenvectors of Ho{X) in S2^\ 
without any spectator-photon 

Figures in the fourth column are the dimensions of the projections of £2^^ on \ i) ® T^"^^ . 

The changes in the levels cq, ei and 62 of Figures 10 and 11 cannot be superimposed 
on one another, in general. Both are to be considered in describing the levels of Hq{X). 
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To complete the study, we have to take an arbitrary number of photons into account. 
The study of the spectrum of Hq{X) is completed when one has also considered eigenvectors 

(2) (2) 

in , , etc.; They give new levels. Eventually, there is a very great number of 
levels. Let S{n,p) be the number of independent symmetric states that can be formed 
with n photons, each photon being in p possible states. Let us consider 1, 2, ■ ■ ■ , n photons 
successively. The fundamental state is split into a doublet, a triplet, • • •, a S'(n, 2)-multiplet. 
If we look for all possible states, all these levels must be considered. As regards level 1, a 
shift, then a doublet, a triplet, ■ ■ ■, a S{n — 1, 2)-multiplet. Lastly, for level 2, we get no 
change, then a shift, a doublet, • ■ ■, a S{n — 1, 2)-multiplet. To be complete, let us recall 
level Co = 0, with eigenvector | 0, O). | 

The doubly infinite degeneracy due to spectator-photons still remains. With the six 
eigenvalues (2,2, ((1,2)1 and ((1,2)2, and lastly ((0,2)1, ((0,2)2 et ((0,2)3, are associated 
eigenvectors 4*^^, ((/)l°^)i and ((/)1°^)2), and lastly ((/)|,°^)i, (4°2)2 et ((/)[,°^)3- As previously, 
through application of (l ® c*{g)^ with g orthogonal to {/oi,/i2}, or more generally 
application of Yl (l®c*((7j))"^", one gets other eigenvectors associated with the same eigen- 
values. In the same way that we built Figure 10 from Figure 9, we could illustrate this 
degeneracy graphically through adding dotted lines in Figure 11. 

Let us come back to our initial problem, which is to determine the spectrum of H{X), 
at least roughly. It may be expected that the perturbation A02 V partially removes the 
degeneracy of each of the six above mentioned eigenvalues, as well as degeneracies of the 
same type, for instance those of the three eigenvalues ((0,1)1, ((0,1)2 and (1^1 that we 
obtained previously. It is this simpler question that we now examine. We are going to 
show that the coupling A02 V splits the first three levels of Figure 10, eigenvalues ((0,1)1, 
((0,1)2 and (1^1 of Proposition 4.2, into an infinity of levels and calculate the splittings of 
((0,1)1 and (1^1, at the lowest order in Ao2- 

IV. 3 Perturbation with respect to Ao2' Second removal of degeneracy 

We are now interested in H{X). The function /02 comes into play. Hence we assume 
that the first vector of the basis (71, (72, ■ ■ ■ of Proposition 4.2 is in Ti^^d- "^^^ fi'i'^' ^ > 1 
are thus orthogonal to H^^^. We saw that the {<i>i1p^'"'^) g = 0, 1, • • •, which are m 
are associated with a unique eigenvalue {Ci,p)k- In the simple cases {i,p) = (1,1) and 
{i,p,k) = (0,1,1), we are going to show that the degeneracy is removed. We calculate 
the approximations of order two in A02 of those eigenvalues of H{X) which tend to ((i,p)fc 
when A02 tends to 0. These eigenvalues depend on q. The approximations are denoted 
by (-2^ p'')fc~^^. The corresponding eigenvectors are denoted by {Xi^p)k and their 2-order 
approximations by (x^ p)^"^'*- We will have 

lim (xS)?^) = (<^S5°'-)), = (1 ® c*ig^)yi4% (4-13) 

(Here, z has not the same meaning as in Section III; // remains zero.) 



Note that, if X ^ Ti,^^^ is an eigenvector of H{X), then for every g orthogonal to Ti 



(3) 
rad' 

^1 <S> {c*{g))^^X is still an eigenvector, associated with the same eigenvalue. This second 
degeneracy removal is thus only very partial. 
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We assume that p = 1 and we limit ourselves to perturbing (Co,i)i) (Co,i)2 if /oi fi2, 
and Ci 1. 

If/oi 7^ fi2, the diagram in Figure 10 transforms into the one of Figure 12. 

^2 (q) |0,1),|0,2>,.. 



'1,1 IN 



Xl,l |1,0>,|1,1>,- 



levels of the 

decoupled 

atom 



levels levels 
of Ho (A) ofH(A) 



(4% 

q=0,l,- 



eigen- 
values 



dimension of the 
field state-space 



q=0,l,--- 



eigenvectors'-^-* and 



(apart from deg.) their content 



Figure 12. Qualitative description of the perturbation of Figure 10 first levels, 
due to the A02 V term; {gi^ = gi) 



The case of (Co, 1)2 is particular as we will see later on. 

Note that the exact eigenvector (xi is no longer in ^^^^ but in ®^o^r+i- 
More precisely, we have 

Proposition 4.4 Through the perturbation A02 V, each of the eigenvalues 
(Co,i)i('^oi) A12) and Cijil'^oij A12) is at least split into an infinity of eigenvalues 
{zq^i)i{X) and -^i^i (A), given at second order in A02 by the following formulas 

(V (4^)P(A) = (Co,i)i(Aoi,Ai2)+A22 (l + n! (^)"' 41i(Aoi,Ai2) (4.14) 
with, if f 01 ^ /12, 



^o?i,i(Aoi, A12) :=n 



A^i |(/02,^l)P 



((Co,i)i)'(A?2(Co,i)i + Agi((Co,i)i-e2)) 



(4.15) 



where |(/o25fi'i)| can be expressed with the Si's, 

and, if f 01 = /12, and thus gi ^ {1 + |soP)"^(/o2 - ^o/oi), 



41i(Aoi,Ai2) 



Aoi 



((Co,i)i)' A?2(Co,i)i + A2,((Co,i)i-e2 



(n|(/o2,^i)|' + ko|') (4.16) 



(a) {zi^i)-^ is obtained through replacing (Co, 1)1 by Ci,i in expressions giving (-Zo'i )l 



(n)x<2 



The residual degeneracy is not mentioned. 
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(in) (Co,i)2j which is more degenerated than the two previous eigenvalues, is also split. 
The 2-order approximations of the perturbed eigenvalues, {zl^{^)f^ , now depending on two 
indices, are obtained through the vanishing of an infinite order determinant. 

The broad lines of the proof are given in Appendix B3, together with the method for 
calculating the corresponding eigenvectors (Xo?i)^^(A), (Xo?i"^'')2 ^('^)' ^^'^ {Xi^i)~'^{^) 

IV. 4 Conclusion of Section IV 

Section IV. 3 described the splitting of the eigenvalues (Co,i)i: (Co, 1)2 and of Propo- 
sition 4.2. It gave a small part of the spectrum of iy(A, 0), described in the second column 
of Figure 12. In view of these results one may reasonably surmise the following points. 

This splitting of the levels of the second column of Figure 9 (or 10) will reproduce 
for those of the second column of Figure 11, which are different. In other terms, the 
degeneracy of the latter levels, due to spectator-photons, will also be removed by the 
coupling X02V. More generally, each level multiplet which was mentioned at the end of 
the proof of Proposition 4.3 is also split when the interaction is totally switched on. This 
rough description of the two level splittings we get by successively taking the two parts of 
the interaction Hamiltonian into account eventually yields quite a complicated spectrum 
for H{X, 0). But the reason of the great number of levels is simple; it is recalled in Section 
V below. 

When fj, increases from 0, we expect that the eigenvalues we found move into eigen- 
values or resonances of H{X,ji). The set of these eigenvalues or resonances thus likely has 
the same rich structure. 

V General conclusion 

Results of Sections III and IV fit into the same frame. They lead us to expect that 
the number of eigenstates or resonances of the atom-photon system is formally, apart from 
accidental degeneracy, the product of the dimensions of the atom state-space and the field 
state-space. The description is complicated because of the multivaluedness of the resolvent 
matrix elements, as functions of the A and /i parameters. 

This picture may be illustrated through the following argument: if the energy of each level 
of the isolated atom is considered as an eigenvalue of -ffatom ® Ifieid, this level is twice 
infinitely degenerated (number of photons and state of each photon). This degeneracy is 
removed with the actual Hamiltonian which describes the coupling of the atom to the field. 
The shift of the naked-atom levels by the coupling of the atom to the photon field is thus 
not the main feature in the change in the Hamiltonian "spectrum". The main feature is 
more the emergence of numerous resonances, as in the monochromatic-photon case. 

For a two-level atom and Hamiltonian (3.1) (Section III), only one photon state comes 

into play. The subspace of Tirad to be considered is 7-^|^^d ^^'^ degeneracies which 

are removed, same energy for states | 0, 0), | 0, ^), • • • , | 0, ^^"■), • • • on the one hand and 
I 1,0), I 1, ^), • • • , I 1, fir^"^), ■ ■ ■ on the other hand, only concern the number of photons. 
A great degeneracy remains since adding photons in states orthogonal to the distinguished 
5f-state do not change the energies. 

For a three-level atom, the number of coupling functions in the Hamiltonian is greater 
and this forced us to start to pay attention to different photon states. As a consequence, 
the just mentioned degeneracy now starts being removed. 
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For a real atom, with its infinity of levels, the splitting will still be greater. Our 

perturbative treatment illustrates how the different photon states may be taken into ac- 
count successively. Calculations will of course be impossible if some physically justified 
simplifications are not made. 

The present limits of the study are the following. 

Section III described resonances for a realistic Hamiltonian 7^ 0), but for a system 
with only two levels. However, even in that simplified case, we are far from having found 
all the resonances since we considered only one-or- two-excitation spaces. Si or 82- It would 
be necessary to take more than two photons into account. But the resonances are then 
given by more and more complicated equations. 

In Section IV, to be able to present a qualitative description of resonances in a three- 
level system, we had to work in the limit fi = 0. This first stage seems unavoidable to us 
if one wants to solve the question completely. Doing this, we were able to take an infinity 
of photons into account. But the calculations are only carried to the second order in A02 
and also the displacement of the resonances when /j, becomes non-zero is just qualitatively 
mentioned. 

However, we have seen that these partial results give new information about hybrid 
states which are present in matter-field interactions such as the interaction to which we 
borrowed our terminology: the interaction of atoms (or molecules) with the electromagnetic 
field. We hope that concrete problems will justify approximations making calculations 
possible. 



Appendix A. Estimation of corrective terms in the Fredholm expansion of 
Proposition 3.4 

Proposition Set (^2(1^, z,p, q) := //^ D2{fJ>, z,p,q) and 





\9ip)\ 


\giq)\ (1 + H)V4 (1 + IqI) 


1/4 (IpI . 






\f{z-i,\p\)\y^ |/(^-/.|g|)|V2 1 


- -21 


\P\\ 


1/2 1 


- -21 


\q\ 


|l/2 1 






\q\ 


12 





\9ip)\y^ 


\9iq)\ 


2/3 (l + |p|)l/4 (l + |g|)l/4 


(N- 


\Q\r 


|/(^-^H)|V3 |/(z-;,|g|)|l/3 


- -21 

M 


\P\ 


|l/3 1 


^ -21 




|l/3 1 


i-\p\- 
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(P3{fx,z,p,q) := 



(P4{f^,z,p,q) := 
defined for z <0. we have 

C2{X, fi, Z) = llf ||(^2(//,^,.,.)l|l (^1) 

ICsiX, fi,z)\< Ms{X, fi, z) := 4 {X^fif \\<p3{fi, z, ., .)| |i {A2) 
Assuming that (3.19) also holds for n = 4, we also have 

\C4X, fi, z)\ < M^iX, fi, z) := 4^ {X^ / ,jif ||^4(//, ^, Olle (^3) 



Proof Use (3.19) and Holder's inequalities. | 
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We omit A, which is fixed to 0.1. Table 2 gives values for C2(//, 2(a*))) 
M3{f^,z^ M), M^{ii, z^ ^il^)) and d;,Ci{ii, z^a\l^)), from 

a,Ci(M,^) = -A^ / 'f^y ^ I 1, (77 — ^—rr^ hn) ^p- 

J f{ii,z- ii\p\){z-2ii\p\)\f{ii,z- ii\p\) z-2ii\p\J 

dp dq 




Pifx, z - fx\p\)iz - 2fx\p\){z - fxip + q)Y 

By definition, Ci (/x, = 1. We see how the terms of the Fredholm expansion 

decrease with the order. 

For n = He -D(/Uc, 0) is close to 0. Let us recall that we have /(A, /Uc(A), O) = 0. It 
would be interesting to see whether I?(/ic) 0) vanishes or not. To try and answer this ques- 
tion, let us estimate the error we made in calculating the zero of D{n, .) by the truncated 
series. Let us consider the following expansion of -D(/x, .) near 2^0, 2(a*) 

1 - Ci{fi,z^ M) - - 4,2if^)) dzCi{fi,zl M) + ^C'2(/i,4,2(/i))+ 

^(^-^0,2(/^)) 5^C'2(/^,4,2(^)) - ^(C3(^,^d,2(^)) + (^-^0,2(^)) ^z^'s 42(A*)) ) 

Let us assume that terms {z — -2o,2(a*))^^C'2(/U, -Zq 2(/^)) ^^id (2 — ZQ 2{f^))dzC3(^fJ,,ZQ 2{fJ')) 
can be neglected. Then the correction to the zero is 

1 



whose principal term is of the order of 4 10 ^. Since it is precisely the order of Zq 2(a*c)5 it 
is not possible to answer the question. 

Appendix B Sketches of proofs of results in the 3-level case 
Bl Sketch of the proof of Proposition 4.2. 

Let $ := I 2,Fn) + \ l,Fn+i) + \ 0, F„+2) be an eigenvector of Ho{X) in £^+2 
associated with z, one of the two eigenvalues (Co,i)i or Ci,i- It can be shown that 

(z - e2)-^A?2(/i2, K+i) V /12 = z-^ XlMu Fn+i) V /oi (Bl.l) 



where 



{fl2,Fn+l){pi,P2,---,Pn) ^ / f 12iP) Fn+l{p, Pi, P2, ' ' ' , Pn) dp 



^ n+1 , ^ 

(r"V/^)bl,p2,•••,Pn+l) = -— tE n. /fe) h{p,) 

1=1 

From (Bl.l), through decomposing F^+i on a basis of the (n + 1) -photon space built 
with /oi,/i2 and the (/^'s, we get that F^+i is a sum of states YliLi with 
Ui — n From relations expressing that $ is an eigenvector, we derive = and, 
if 2 ^ 0, F„+2 = z~^^Jn + 2 Aoi-Fn+i V /oi; this implies that $ is in ^0,1,1 if 2; = (Co,i)i 
in Qi^i if 2; = Hence (i) and (iii). 

If z = (Co,i)2 = 0, then $ is an eigenvector if and only if F„ = 0, = and 

(/oi,i^n+2) = 0. Hence (h). 
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B2 The six eigenvalues (2,2, ((1,2)1, ((1,2)2, ((0,2)1, ((0,2)2 and ((0,2)3, with 
eigenvectors in S^^ . 

The space spanned by vectors 

I 2,0), |l,/i2), |l,/oi), |0,/i2V/oi), |0,/oiV/oi) 

is invariant. Through adding | 0, /12 V /12), we get a basis of £2'\ in which the matrix of 
Hq is 



/e2 


A12 


S1A12 











\ 


A12 


ei 





V2 





\/2siAoi 










ei 




V2A01 










\/2Aoi 






















v^Aoi 













V 




















One of the eigenvalues is and the eigenvector is the one mentioned in the text. The other 
eigenvalues are the ('s for which 



((( - ei) (((( - ei)(( - 62) - 3A2,(C - 62) - A^^C) + C>^l^{2Xl^ + (2 - \s^\'')Xl^) - 2e^X 



4 
01 



vanishes. Neglecting fourth order terms, we get the following solutions at second order in 
A 

((0,2)2 = 0, ((0,2)3 = Aqi, ((°2)i = ei, (("2)2 = ei H Aqi H — A12, 



ei - 62 



(2° 2 = 62 + 



-A 



62 - ei 
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This proves that these eigenvalues are different from those found in 8.^^ . 
B3 Sketch of the proof of Proposition 4.4. 

Let us use Kato's method^"^ to determine the three (infinite) sets of perturbed eigen- 
values which tend to each of the unperturbed eigenvalues. (See a short account in Ref. 
24.) 

(i) Let us first consider ((o,i)i- The unperturbed eigenspace is ^0,1,1 (Proposition 
4.2). Let Pq'^'^ be the projector on this space and Qq'^'^ = 1 — Pq'^'^. We need operator 
Qo'^'^ which is sometimes written as Qo'^'''^[((o,i)i — -f^o]~^Qo^'^; for x E it is defined 



0,1,1 



X. 



by Qq ' X = Qq' ' z, where z is any vector in S satisfying [(Co, 1)1 ~ Hq] z = Qq 
Let '^^o,l,l be the direct sum of the eigenspaces associated with eigenvalues of ^^(A) which 
tend to (Co, 1)1 when A02 tends to 0. Let P^'^'^{Xo2) be the projector on 7io,i,i- There- 
fore P°'i'HAo2) — > pS'^'^. We assume that Pq'^'^ and P°'"'^'"'^(Ao2) establish one-to-one 

-^02^0 

correspondences 7Yo,i,i — ^ ^0,1,1 and ^0,1,1 ^ ■^0,1,1- Set 
L 
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0,1,1 = pS''''HP'''''pS'''' and Ko,i,i = pO'MpO,i,ipO,i,i 



(B3.1) 



We recall24 the expansion of P^'^'^: with ^(o) := -P°'^'^ and S^''^ := (Qo'^'^)'=, 



pO,i,i ^ pO,i,i _ ^ AJJa 5] 5(feo)y5(fei) . . . (B3_2) 

n A;i>0, feo4 l-fe„=n 

A necessary and sufficient condition for x £ ^0,1,1 to be an eigenvector of H{X) associated 
with the eigenvalue z is that there exists G ^0,1,1 satisfying x = P^'^'^4> and 



-0,1,1 



= z Ko,i,i <p (B3.3) 



-f'0,1,1 and Kq i i are operators in ^0,1,1 and the problem of the perturbation of (Co, 1)1 is 
turned into finding such z. We still have a degeneracy due to photon-states in ('^rfd)^' 
indeed, if is a solution for (B3.3), then ni>2(l(>?)c*((7j))"'"(^ is still a solution. Thus we are 
not going to look for all ^'s, but only for those in an invariant subspace of ^0,1,1 • Lemmas 
B3.1 to B3.6 prepare the calculation of eigenvalues into which (Co, 1)1 splits. The result for 
Ci,i will be obtained through a simple change in the notations. The splitting of (Co, 1)2 is 
just outlined. 

Let gi be the function of hI^J^ orthogonal to Ti^ad ^0^1, 1 subspace of 

^0,1,1 spanned by (l <S)C*{gi))^ (<?^o°i)i5 > 0. Let us denote the approximations of -f^o,i,i 
and 1/0,1,1 to order q in A02 by K^J -^ and -£'0,1,1- 

Lemma B3.1 For all q, Qo]i^i is invariant by K^J -^ and -£'0,1,1- 
Proof Operators 1 <S> c*{gi) and V sent S^^^ into since functions orthogonal to "^^^ad 
do not play any part. Now, Pq'"^'"^ sends £^^^ into ^o^i^i ^^^^ ■ Thus, if X is any 
endomorphism of £^^\ '^'^XPg '^'^ is an endomorphism of Qq^I Qq'^'"^ also leaves S^^^ 
invariant. The same is true for Qq'^'^ and therefore for P^'^'^. Hence the Lemma. | 

Expressions to the second order in A02 of -Z-0,1,1 and -f^o,i,i are 

^o-M = pS'''' - >^l2PS'''' V (g^'^)' V pS'''' (B3.4) 
^o-li = (Co,i)ii^P + Xl, pS'''' V QS'''' V pS'''' (B3.5) 

To calculate -£^(^i,i(0o?n-i-i)i -^0,1,1 (^o?n+i)i' need the following lemma: 
Lemma B3.2 

^(C+i)i= (B3.6) 
with {(fi^o})i = Ni j^^so andjorn > 1, 

(^'S)! = {so gr+n{fo2,gi) /oi V^r") (B3.7) 

[^0,1)1 



w;/iereAri = (l + (Co,i)r'A§i)-V2 
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Lemma B3.3 K^J -^ and L^^ i are diagonal in the basis (0o"n+i)i ^/ ^0^1,1 • 
Proof ^((/>o^^_|_i)i G since (0o^^+i)i has no component on | 2) ® Through using 

<5S'^'^^n+2 C ^^+2: we then get 

and ii'Q-;^ ^ is diagonal. The same is true iov L^^ ^. | 

Lemma B3.4 Let us set ui^n '■= Qo'^'^ V (^o"n+i)i ^'^^ define Ki^n, ^i,n+i <ind hi,n+2 
by wi,n = I 2, Ki^n) + I 1, 6*1,^+1) + I 0, /ti,n+2)- We denote the expression of z to the 
second order in A02 by {zQ^i)f'^ . We have 

(4^ )P = (Co,i)i + A^2 «l,n) (B3.8) 

Proof From Lemma B3.3, it follows that eigenvectors of i and K^i 1, 

associated with eigenvalues which we denote by (/o"'')P (^o"^^)P respectively. Vectors 
(f) in ^0^1 1 satisfying (-^0,1,1 "~ ^^0,1,1)'^ — necessarily these ((/'o^i+i)i; the corre- 
sponding 2-value, for each n, is ((4"i)p) ^(4'\'')P- Through using (B3.4) and (B3.6), 
we get 

(45 )P = ((<^£+i)i,^o-?,i(4i+i)i) = 1 - A^2ll«l,n|P (B3.9) 

aS)P = (Co,l)l(l - Ag2||«l,nir) + Ag2(V^ 

As a consequence, 

(4?)P = (C0,l)l(l - A^2ll«l,nir) + A^2((<^S)l,«l,n) (B3.10) 

Hence (B3.8) holds since 1 — | |wi n| P is to be replaced by 1, to the considered approximation. 
I 

Lemma B3.5 

(a) Let us set Mi :— ( ^"-ea Tcfih) ' n > 0, a vector v satisfying 

[(C0,l)l-^0]t^= |2,(^g)l) ZS^;i,n= I 2, \^,K^n+2),Wlth 



«in = 77 — ^ hVn+lAi27- — 7 (B3.11) 

' (Co,i)i-e2 (Co,i)i-e2 

/ii,n+2 = Vn + 2 Aoi — 7—^; (B3.12) 

^i,n+l = iVlx/^ ..^ ^ (^O(^;,n+l)^^^+n(/o2,ffl)(^;,n+l)^^0 (B3.13) 

((Co,i)i -e2)(Co,i)i 
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where, z//oi ^ /12, 

{o[,n^.r' = ^^"■^]r'% I |. (^1/01 - /12) v,r (b3.i4) 

^^12 ~ |-^l| 

(B3.15) 

If foi = /l2j 

(^'i,n+i)^'^ = -^1 /oiV^r (B3.16) 
(^m)^'^ = 0' i(^i,n+ir^ - -\m, foi V /oi V (B3.17) 

r&y' ^1,71 of Lemma B3.4 zs equal to vi^n 
Proof Relations (B3.11) et (B3.12) are obtained through projecting the equahty defining 
V on I 2) ^Hrad and | 0) <8)Hrad) respectively. The projection on | 1) (8)Hrad implies that 
^i,n+i satisfies 

A,n+l^'l,n+l = V/12 (B3.18) 

where 

'Cl,n+l^i,n+l '■= -77 r^-^f 12 V (/l2,6'i,n+l) " ^7 — foi V (/oi, ^'l,n+l) (B3.19) 

Uo,iji-e2 {^0,1)1 

For n > 1, (<^o?i)i, given by (B3.7) is decomposed into two parts. Hence we introduce two 
functions (d'l^n+i)^^^ and (^i,n+i)*'^^ satisfying 

A,n+l(^i,n+l)^'^=^rV/i2 (B3.20) 

A,n+i(^l,n+i)^'^ = gi^""'^ V /oi V /12 (B3.21) 

so that a solution of (B3.18) will be given by (B3.13). To prove (a) if /oi 7^ /12, we check 
that (B3.14) satisfies (B3.20) and that (B3.15) satisfies (B3.21). We proceed in the same 
way with (B3.16) and (B3.17), if /oi = /12, (B3.18) and (B3.19) being stiU true. Only 
(^i.n+i)*'''^'' plays a part if n = 0. 

To prove (b), we note that vectors in ^0,1,1 are linear combinations of | l,gi-^ V 
• • ■ Qi^^) and I 0, foi y Qi^ V ■ ■ ■ (7i„). Since (a) implies that O'l is a sum of symmetric 

/OA 

products of terms one of which at least is in 'H).^^, we have {O'^ ^j^i.gi^ V ■ ■ ■ Qi^^^) — and 

(^'l,n+2: /oi Vfifii V ■ ■ 

■ 9in+i) 0- a consequence, | 2, Hi^^ 

,n+i) and I 0, h'i^^_^2) 

are orthogonal to ^0,1,1 and ui^^ = (5o'^'^fi,n = vi,n- I 
Lemma B3.6 Take n>0. For /oi ^ /12, 

fr,An)\ ^ \ ^Ar2 Aqi |(/o2,fi'l)P rvtrt 00\ 

{m,i)u «i,n) = nN^ -— — "Y (B3.22) 

((Co,i)i)'(A?2(Co,i)i + A2,((Co,i)i-e2)) 
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For /oi = /i2, 

((<^S)i,«i,n) = iVf- -—- ^t^— ^(n|(/02,^i)|^+|.0|') 

((Co,i)i)'(A?2(Co,i)i+A2i((Co,i)i-e2)) 

(B3.22') 

Proof Use (B3.7) and Lemma B3.5. | 

Formulas (4.14) and (4.15) then foUow from (B3.10) and (B3.22). In the same way, 
(4.16) follows from (B3.22'). 

(ii) All that has been written up to now can be transposed from (Co,i)i to with 
the following changes 

Ni, Ml, (<^o"i)i are thus also changed, as well as £i,n+i and the solutions ^i,n+i 

/ll,n+2 and Kl^n- 

One then get (ii) of Proposition 4.4. 

(iii) Let us come on now to the splitting of (Co,i)2) an eigenvalue which is zero and 
exists only if /oi ^ fi2- The issue is more complicated due to the fact that the unperturbed 
eigenspace of interest is no longer just spanned by the 11(1 ^ c* (gi))"^ {(/) 0*1)2 ■ It is spanned 

/ \ ^2 12 

by the | 0,n(c*((7i)) 'O), with i = or 1. Let us denote this space by Q. Let -Pq ^® 
the projector on With ^ ^-i^d L'^\ ^ defined as in (B3.4) and (B3.5), we can see that 

2 has non- vanishing matrix elements between say g^^ V (7^ and 5^9 ''^"'"^^ 9\^^^'^^ 1 
or between \J (^l^^^^^ and ^"^^"^^ g^x^^^ °^ between \J g^^^^^ ^^"^ 

^v(p 1)^ This makes the computation of the satisfying det(L^^ 2 ~ -^^ ^012) — more 
intricate and we don't calculate them here. 

This completes the proof of Proposition 4.4. | 

Eigenvectors (Xo"i)i^\ (Xi"i)^^'' obtained from the correspondences 

^0,1,1 ^0,1, 17 and Qi^i —>■ Ti.i,i, through the second order expansion of operators 
P'^'^'^(Ao2) and P"'^'"'^(Ao2) which perform these correspondences. For example, we get (see 
Ref 24 p. 614) 

(xS)f ^ - + A^2 ( - Po°'^' V(Qr '^)^^ + Q^' Vq^' V) (</>£+i)i 
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